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I. INTRODUCTION 

The understanding of the structure of bound state 
systems in nuclear and particle physics requires the 
development of a relativistic nonperturbative frame- 
work. For obvious practical reasons, any calculation 
of this type relies on approximations, and one thus 
needs a systematic strategy in order to improve the 
approximations which are made, in complete analogy 
with perturbation theory. 

Light-front dynamics (LFD) is a very powerful tool 
to calculate bound state properties. It is one of the 
three forms of dynamics proposed in 1949 by Dirac [T]. 
In order to perform systematic calculations of physical 
observables on a large scale, one should however be 
able to solve three important problems. 

The first one is the explicit violation of rotational 
invariance by the choice of a given light-front plane. 
The control of this violation is important in order to 
define unambiguously all physical observables. This 
can be done simply using the covariant formulation of 
light-front dynamics (CLFD) [IE]. 

The second one is the possible appearance of un- 
canceled divergences when the Fock space is trun- 
cated, in any approximate nonperturbative calcula- 
tion. One thus should make sure that no divergences 
are left uncanceled. This is enforced using the Fock 
sector dependent renormalization scheme |4j. 

Finally, one should develop a regulariza- 
tion/renormalization scheme which preserves all 



symmetries, and which is well adapted to extended 
numerical calculations. 

We have advocated in a previous study [5] the use of 
the Taylor-Lagrange renormalization scheme (TLRS) 
0. We have shown in particular that this scheme 
is very well adapted to any calculation in LFD. It 
is systematic, can treat singularities of any type on 
the same footing, and moreover does not require to 
perform any infinite scale limit. 

Many other regularization methods are available in 
the literature. Let us mention here the most impor- 
tant ones. 

i) The cut-off method. This is a simple, but to some 
extent brutal, way to regularize divergent amplitudes. 
It however violates gauge invariance and, in LFD, ro- 
tational invariance [7]. It should thus be avoided in 
any realistic calculation. 

ii) The Pauli-Villars (PV) subtraction method. In 
LFD, the PV regularization scheme amounts to ex- 
tend the Fock space to include PV fields with negative 
norm [1]. While this method is attractive for its sim- 
plicity and immediate use, it has some serious disad- 
vantages in any systematic calculations. The number 
of PV fields may be large if singularities are of high 
order, as in effective field theories. This may imply a 
large number of PV components which are not easy to 
implement in systematic calculations. Moreover, the 
limit of large PV masses has to be performed numer- 
ically. This may also be rather delicate to achieve in 
large systematic numerical calculations. 

Hi) The dimensional regularization method. This 
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method is largely used in perturbation theory in the 
covariant Feynman approach. In LFD however, it has 
never been used since one would have to reformulate 
LFD in arbitrary D dimensions. This is also not in the 
spirit of LFD which deals only with physical degrees 
of freedom. 

iiv) The Bogoliubov-Parasiuk-Hepp-Zimmermann 
(BPHZ) method. In this method each Feynman graph 
is treated separately. The contribution is made finite 
by subtracting as many terms as necessary from the 
Taylor expansion at zero external momenta of the 
integrand. All Feynman integrals being convergent, 
no intermediate regularization is required, thereby 
showing that all n-point functions are effectively 
regularization independent. They only depend on 
the renormalisation conditions. However there are 
some complications for zero masses and despite its 
theoretical importance, the BPHZ scheme is not 
easy to deal with in practical calculations. The link 
between the TLRS and BPHZ schemes resides in the 
use of specific test functions equal to their Taylor 
remainder of any order. 

We shall exhibit in this study the relationship be- 
tween TLRS and the PV-type regularization proce- 
dure. We shall illustrate this relationship in the sim- 
ple, but instructive, calculation of the various com- 
ponents of the self-energy of a fermion in second or- 
der perturbation theory in the Yukawa model, paying 
particular attention to the restoration of rotational 
invariance. 

The plan of the article is the following. We recall in 
Secinithe main features of the TLRS rcnormalization 
scheme. We then calculate in Sec. |III| the self-energy 
in second order perturbation theory in the Yukawa 
model. We draw our conclusions in Sec. HVl 



II. THE TAYLOR-LAGRANGE 
RENORMALIZATION SCHEME AND 
PAULI-VILLARS SUBTRACTION 

It is a common lore fH* that any field (taken 
here as a scalar field for simplicity) should be consid- 
ered as an operator-valued distribution. This means 
that it should be defined by its application on test 
functions, denoted by p, with well identified mathe- 
matical properties. In flat space, the physical field 
if(x) is thus given by [S] 



in D dimensions. If we denote by / the Fourier trans- 
form of the test function, we can further write ^p{x) in 
terms of creation and destruction operators, leading 
to 



(27r)(D-i) 2£p 



(2) 

with = + TO^ . 

From this decomposition, it is apparent that test 
functions should be attached to each fermion and bo- 
son fleld. Each propagator being the contraction of 
two fields should be proportional to p. In order to 
have a dimensionless argument for /, we shall intro- 
duce an arbitrary scale A to "measure" all momenta. 
A can be any of the masses of the constituents. To 
deal with massless theories, we shall take some arbi- 
trary value. The final expression of any amplitude 
should be independent of A. In CLFD, the test func- 
tion is thus a function of ^ only. 

As recalled in 0, the test function / should have 
two important properties: 

i) The physical fleld f{x) should be independent of 
the choice of the test function. This later should there- 
fore be chosen among the partitions of unity (PU) . It 
is a function of finite support which is 1 everywhere 
except at the boundaries. This choice is also neces- 
sary in order to satisfy Poincare invariance since, if 
/ is a PU, any power of /, /", is also a PU. In the 
limit where the test function goes to 1 over the whole 
space, we then have /" — > / and Poincare invariance 
is recovered. 

ii) In order to be able to treat in a generic way 
singular distributions of any type, the test function 
is chosen as a super regular test function (SRTF). It 
is a function of finite extension - or finite support - 
vanishing with all its derivatives at its boundaries, 
either in the ultraviolet (UV) or infrared (IR) domain. 

Any physical amplitude is thus written in a 
schematic way like 



A 



dX T{X) fix) 



(3) 



ip{x) 



d y(j){y)p{x-y) 



(1) 



for a one-dimensional distribution. In this form, the 
amplitude A does not differ from the calculation using 
a cut-off procedure. In the UV domain for example, 
the cut-off, denoted by H , would correspond to the 
support of /, with f{X > H) = 0. In order to go 
beyond the use of a naive cut-off, we should investigate 
the scaling properties inherent to the limit X ^ oo 
since in this limit rj'^X goes also to oo, where 77^ is 
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an arbitrary dimensionless scale. To do that, we shall 
consider as an example a distribution T{X) leading 
to an X integral diverging like log(iJ) in the absence 
of f{X) and use the Lagrange formula written in the 
following form, in the UV domain : 



with 



f{aX) - -X 



dt 
T 



dx [f{Xt)] , 



(4) 



for an arbitrary intrinsic scale a which can be chosen 
positive if T has no other singularity at finite X . 
This formula is an identity for any function / which 
is a SRTF. In order to introduce the arbitrary scale 
77^, we shall consider a running boundary condition, 
i.e. a boundary condition which depends on the given 
variable X using 



H{X) EE fj'XgiX) 



(5) 



up to an additive arbitrary finite constant irrelevant 
in the UV domain. Note that the support of the test 
function is the same in the right- and left-hand sides 
of Eq. 0. This implies that Xt < H{X) for any 
argument Xt of the test function (see [5] for more 
details) . 

The variable X on which the running condition ([5| 
is applied should not be linked to any intrinsic scale. 
This is necessary in order to make sure that the limit 
X — ^ 00 is properly done, i.e. that X should be larger 
than any other physical scale present in the amplitude. 
Once this is done, the test function depends generally 
on aX, where a is a priori a function of the kinematical 
variables of the system under consideration, and one 
should consider the Lagrange formula in the form Q . 

In order to extend the test function to 1 over the 
whole space, we shall consider a set of function g{X), 
denoted by ga{X), where by construction a is a real 
positive number less than 1. A typical example of 
ga {X) is given in [S] , where it is shown that in the limit 
a — ^ 1~, with rf > 1, the running support of the PU 
test function stretches then over the whole integration 
domain, and / — >■ 1. In this limit ga{X) — > 1. 

After integration by parts, the amplitude A writes 



I dX dx [XT{X)] 
la 



t 



f{Xt). (6) 



In the limit a — !■ 1 , the requirements are such that 
ga{X) — 1 and f ^ 1. One can thus define the 
extension in the UV domain, denoted by T^, of the 
singular distribution T by 



A 



dX T>{X) 



(7) 



T>iX)=dx[XTiX)]Log 



r 21 

7f 



(8) 



where the derivative should be understood in the sense 
of distributions [5J . It depends logarithmically on the 
arbitrary scale 77^, with 77^ > 1. The amplitude ([t]) 
is now completely finite. Note that we do not need 
the explicit form of the test function in the derivation 
of the extended distribution T>(X). We only rely 
on its mathematical properties and on the running 
construction of the boundary conditions. 

The running boundaries of the test functions are es- 
sential for the preservation of symmetries which would 
otherwise be destroyed with the usual cut-off test 
functions. Qualitatively, the reason for this property 
is the following: the cut-off functions are equal to 1 up 
to a point Xc and fall then down to zero over a finite 
interval with some shape which does not change when 
Xc is sent to infinity. On the contrary, with running 
boundaries, the width of the region where the test 
function falls from 1 to zero increases proportionally 
to Xc when Xc goes to infinity, implying an infinites- 
imal drop-off of the test functions in the asymptotic 
limit. We call this behavior an ultrasoft cut-off. 

The extension of singular distributions in the IR 
domain can be done similarly [S]. For an homoge- 
neous distribution in one dimension, with T[X/t] — 
t^^^T(X\ the extension of the distribution T in the 
IR domain writes ^ 



T<(X) = (-l)''c)i+^ 



k\ 



T{X)hog{f,X) 



(9) 



The extension [X] differs from the original distri- 
bution T{X) only at the singularity ai X = 0. 

The amplitude ^ can also be transformed alter- 
natively in order to exhibit a PV-type subtraction. 
Using the Lagrange formula in the form 



f[aX] 



dt dtf [Xt] 



(10) 



^ Here and in Eq.|9]l below the extensions for T{X) are only 
valid when taken under the X-integral symbol, see [5] for 
discussions 

2 At d dimension, k is defined by T[X/t] = t^+'^T{X). This 
corrects a misprint in [5]. 
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we can rewrite the physical ampHtude A in the fol- 
lowing form, after the change of variable Z — Xt and 
in the limit a — )■ 1^ 









/ dZ dt dt 






Jo J a 







A 



With a typical distribution T{X) = 
intrinsic scale a, one thus gets immediately 



(11) 



with an 



T>{X) 



1 



X 



X 



77^ 



(12) 



We recover here a PV-type subtraction, with a scale 
7]^ . This scale is completely arbitrary, with 77^ > 1 and 
not compulsory infinitely large as required for the PV 
masses. This extension of T{X) leads to a well defined 
amplitude Q. 



III. APPLICATION TO THE CALCULATION 
OF THE SELF-ENERGY IN LIGHT-FRONT 
DYNAMICS 

CLFD was first proposed in Ref. and detailed in 
the case of few-body systems in Ref. |3j . Any physical 
system of momentum p is described in LFD by a state 
vector projected onto the plan = t + ^. In CLFD, 
the state vector is defined on the general plane deter- 
mined by the equation a = lo-x, with = 0. The 
covariance of our approach is due to the invariance 
of the light-front plane equation. This implies that 
ui is not the same in any reference frame, but varies 
according to Lorentz transformations, like the coordi- 
nate X. It is not the case in the standard formulation 
of LFD where to is fixed to cj = (1,0,0,-1) in any 
reference frame. The evolution of the system is thus 
defined in terms of the light-front time a. 

We shall consider in the following the simple case 
of the self-energy of a fermion in the Yukawa model, 
in second order perturbation theory. From a practical 
point of view, any amplitude is calculated using the 
equivalence of Feynman rules, as detailed in Ref. [5]. 

In second order of perturbation theory, the self- 
energy S(p) (up to a conventional minus sign) is de- 
termined by the sum of the two diagrams shown in 
Fig-E 



E(p) = T,2bip) + ^fc{p) 



(13) 



They correspond to the two-body contribution and 
the fermion contact term, respectively. These dia- 
grams correspond to time ordered diagrams in the 



light-front time a. Analytical expressions for the cor- 
responding amplitudes read 



5^26 (P) 



9 



S'^^\p + CUT - kl - k2) 



dT 



(27r)3 

{jii + m)e{uj-ki)S{kl - m^)f 



A2 

2 



A2 



d^ki 



d^k2 , (14) 



^fc{p) 



9 
(2n) 



{u:-k2)5{kl - fi^) 
2uj-{p - k2) 



A2 



(15) 

where g is the coupling constant of the fermion-boson 
interaction, and m and /j, are the fermion and boson 
masses, respectively. Since the test functions are PU, 
we shall identify in the following with /. All the 
particles are on their mass shell in LFD, but off en- 
ergy shell, so that the momentum r represents the off- 
shell energy of the intermediate fermion-boson state 
The momentum p corresponds to p ~ pi ~ loti 
where ti is proportional to the off-shell energy of the 
initial, or final, fermion, and pi is the on-shell four- 
momentum of the fermion. 



(a) 



FIG. 1: Two contributions to the LFD fermion self-energy 
— E(p): the two-body intermediate state (a) and the con- 
tact term (b). The solid and dashed lines represent the 
fermion and the boson, respectively. 



Note that in LFD, the self-energy may depend a 
priori on the position of the light-front characterized 
by Lj. In an exact calculation, or in perturbation the- 
ory, we should check explicitly that this dependence 
disappears in order to recover the well known results 
in the 4D Feynman approach. 

The self-energy can therefore be decomposed in the 
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most general spin structures according to 



Mi/i 



M 



ui-p 



(16a) 
(16b) 



where the coefficients B, and C are scalar functions 
which depend on only. They are inde pen dent of 
uj. We denote by M^. The scale M in ( |16| is just 
introduced for convenience in order to have the same 
dimension for all the coefficients A,B,C and C/c- 
These coefficients can easily be calculated according 

to 



g'A{M^) = ^Tr[I](p)] , 
g^B{M^) = 



M 
1 



(17a) 

Tr[E(p) ^] , (17b) 



AM' 



-Tr 



S(p) ^ 



LJ-p 



.(17c) 



The coefficient Cfc is a constant. It can be extracted 



directly from Eq. ( 15 ) 



In order to transform 
the Lagrange formula 



14| ) using TLRS, we shall use 
4 ) in a slightly different form 



POO 

f [a^X^] =- dt dtf [XV] , 

J a 



(18) 



where the scale a should be identified later on. With 
the change of variable fci = kist, k2 = k2st, and f = 
Tst, we have 



^2b{p) = - 



(27r)3 



d*ki 



df 
f 



dt dt 



ds 9s— S'^'^^pst + ojf — fci — ^2) 



1^2 



e{uj-k2)Kkl - M^st)(/^i + mst) e(oj-ki)5{kl - m^st) f 



2A2 



t2a2A2 



(19) 



The four-momentum conservation law, and the on- 
mass shell conditions in Eq. ( 19 ) are equivalent to the 

(20) 



original ones in ( 14 ) after the transformation 



m, /i, p — >■ mst, fist, pst . 
Using the kinematical variables defined by 

UJ-k2 



ui- pst 



, R = k2~xpst with R = (i?°, R_l, i?" 



(21) 

we have, in the reference frame where p = 0, R_l = 
k_L. Since uj-R = 0, we also have i?° = i?". The 
momentum kj^ is the perpendicular component of the 
four-momentum k2 with respect to the position uj of 
the light-front. With the transformations (20), we 
have thus, in the limit of large momenta (UV regime) 



and using the kinematics detailed in Appendix B.l of 

151 



k2 



kl 



4x2M2s2t2 

kl 



4(1 - a;)2M2s2t2 



(22) 
(23) 



To simplify the notation, we shall define the dimen- 
sional scale 



a{x) 



ni^x + fi^{l-x)- APx{1 - x) 



(24) 



With the change (20 1, the coefficients A,B and C 
thus write [7] 
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167r2 



Mg2 
'l67r2 



J a 



QQ pOO 



dx / dt df ds d 



+ 'm?sH'^a{x) 



-'a 



CXD fOO 



dk]^ / da; / dt dt j ds ds r-^ 



(l-x) 



k^ + m'^s'^t'^a(x) 



/[]/[], 



/[]/[] , 



J a 



OO /'OO 



dki / da; / dt 9* / ds d, 



52^2 (1 - x) [k5_ + m'^s^t'^a{x)] 



/[]/[ 



(25a) 



(25b) 



(25c) 



where the notation /[]/[] stands for 



/[]/[] = / 



4(1 



2M2a2A2s2t4 



(26) 



The additional factors st and in (25b) and (25c 



respectively, as compared to (19), originate from the 



momentum dependence of the self-energy in the de- 
composition ( 16a) with the replacement p — ^ pst. 



In order to calculate the coefhcient A, we can sepa- 
rate the integration over x in two domains, for x < ^ 
and X > |, as already done in [5]. In the first domain, 
we can first eliminate the integral in s by r ede fining 
k_L = — , and using the Lagrange formula (18 1 on s 
backward. It thus reads 



dt dt 



167r2 



dkj^ / da; 
Jo 



k^ -|- iv?t'^a{x) 



kl 



4x2M2a2A2t2 



(27) 



where the second test function has been put to 1 since 
its argument is always smaller than the one retained 

in (27). With the change of variable Z ~ 2xMAt ' 
get 



By eliminating the intrinsic scale a (a;) with the iden- 
tification a = a{x), and with Z = aY , we get 



A{M^) 



167r2 



dY 

"'0 

2a;MA 



da; 



dt dt 



a.{x) 



m? -i^Y + t 



f[Y^] . (29) 



The integration over t gives simply, with / — > 1 and 
the upper limit fixed by the running condition y2^2 ^ 
H{Y'^), i.e. t < rj 



A{M^) 



mg 
16^ 



dY 



2xMA 





1 

23:MA y 



da; 



2xMA 



Y + a{x) 



(30) 



A{M^) = 



2xMA 1 

^2 + 



dt dt 



Z2 



The calculation in the interval ^ < a; < 1 and a first 
(28) integration over t instead of s gives the same inte- 
grand, so that we finally get 



7 



,2 



Jo 



mq , mq 

■ — — log 7] H — 

167r2 ^ ' ^ 16^2 



1 



1 



2xMK 



Y 



dx log 



-7? 



2 a; MA 



y + Q:(a;) 
- APx{1 - 



(31) 



This result is the same as the one already given in 
[5] . The calculation of B proceeds in exactly the same 
manner, with just an extra factor (1 — x) in the inte- 
grand. 



We shall now concentrate on the calculation of the 



coefScient C. Using (25c I, we can decompose C in two 
parts 



C(Af2) 



32tt^M 



9 



327r2M 



dkl 



dkl 



dx / dt df ds d. 



m?dxOi{x) 



+ Tm?s'^t'^a{x) 



/[]/[] 



dx 



dt dt 



(32) 



C1+C2 



Let us first calculate Ci. With the change of variable 
a{x)s'^t'^ — u, we have 



Ci = - 



9 



32n^M 



00 /"OO 

dt dt / ds 9, 



1 



dki 



m?du 



M^.^t^ k2 

2 -L 



-/[]/[]■ (33) 



The test functions provide the convergence of the in- 
tegral in k^, as well as an upper limit, 77, in the t and 
s integrations from the running condition on the test 
functions. The order of integrations can be changed 
at will and we remark that it is legitimate at this stage 
to set the test functions to 1 since, after integration 
over u the integral in k^ is henceforth finite. We thus 
get 



327r2Af 



dt dt I ds d. 



dk\ 



log 



log 



sH^ 



(34) 



After a change of variable X = — j + 52^2 ^^j^g gj-g^- 



term, and X 



s t in the second, we get 



? 2 



dt dt I ds ds 



327r2M ./„ \l„ " sH^ 
which finally gives, with X = YsH"^ 



dX log X . 
(35) 



2 2 
g'^m, 

'32^ 



2 rn i-ri 1-1 

— j dtdtjds ds j^^ dY\og[Ysh''] 



(36) 

The final integrations over s and t lead to a double 
PV-type subtraction 



2 2 
g'^m 

327r2M 



log 



dY 



a2 



log[r] - log 
+ log 



Y 



411 

a4 



. (37) 



This result is very similar to the calculation of the 
coefficient C in [7] where it was shown that it is zero 
with one PV fermion and one PV boson subtraction. 
These subtractions are here provided by the integra- 
tion over the variables s and t. Note that we indeed 
need both subtractions to get the final result. 
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This result relies on the property that the numera- 
tor of the integrand in Ci in (32) is just the derivative 



of a{x). This is not the case for any other compo- 
nents. Because of this peculiarity, the calculation of 
Ci proposed in [5] (coefhcient J2 in Appendix B.3) 
is only correct for this case since otherwise it would 
lead to a null result independently of the form of the 
singular distribution. In the calculation of 5 , the 
identification of the product of the two test functions 
with only one symmetric in the change x ^ 1 — x is 
legitimate since under this change only the boundaries 
of the M-integral are interchanged thereby giving zero 
for J2. 

We turn now to the calculation of the coefficient C2 ■ 
Since the integrand has no intrinsic scale, it is natural 
to calculate C2 using the IR extension given in ^ 
for an homogeneous distribution. With the change of 

variable — y , we can integrate over s and t using 
the Lagrange formula Q back. We get 



C2(Af2) = 



dY 



32tt^M Jo Jo 1-x 



dx 



f 



1 



43.2^271^2^2 



/ 



1 



4(1 - x)^Y^]VPA^ 



(38) 



To keep the symmetry a; — )• (1 — a;) in the integrand, 
we shall rewrite C2 as 



C2(M2) 
/ 



327r2M 



1 



42-2^2^/2^2 



/ 



°° dY n- dx 

Y^ Jo x{l - x) 
1 

4(1 - 2;)2y27v/2A2 



(39) 



In the domain x < 1/2, the first test function only 
matters, i.e. is different from 1 within the relevant 
boundary. We have therefore, with the change of vari- 
able 2YxMA = u 



C2(M2) 



167r2 



dx 



du 



1 



(40) 



The test function just provides the extension of the 
distribution at m = from ^ , which is the pseudo- 
function of 1/^2 [5]. We finally have 



C2(M2 



log 2 



16^2 



du Pf 



0. (41) 



We find that the coefficient C{M'^) is identically zero, 
as it should be. 



The calculation of the contribution from the contact 
interaction can be done very similarly to the calcula- 
tion of C2. We have, from (15) and (16b I [7 



C 



fc 



9 



327r2M 



dki 



dx 



x{l — x) 



f 



Ax^M^sH^ 
(42) 

The integration over x can be decomposed in two 
parts, for x < ^ and for x > ^. In the first one, 
with = ^, we get Cfc = -C2{M^) = 0. The 
second one is identically zero with a principal value 
prescription to calculate the integral at x=l. This in- 
sures that the self-energy S(p) is indeed independent 
of the arbitrary position of the light-front. 



IV. CONCLUSIONS 

We have shown in this study that the recently pro- 
posed TLRS [5] leads naturally to a regularization of 
any amplitude very similar to a PV-type subtraction. 
However, contrary to the original PV regularization 
method, the TLRS does not necessitate to perform 
any infinite mass scale limit. 

The application of our formalism to the calculation 
of the self-energy in second order perturbation theory 
in the Yukawa model, using LFD, is very instructive. 
It is very similar to the standard calculation using 
PV fields with a negative norm 7|. The coefficients 



A and B in the spin decomposition (16a) depends in 
TLRS on an arbitrarily dimensionless scale 77, while 

it depends on in the PV method, where Apv is 
the PV boson mass. This mass must be taken very 
large compared to any physical mass scale present in 
the amplitude. In both methods, the regularization 
of the amplitude is achieved by using a single PV- 
type subtraction. On the other hand, the coefficient 
C is identically zero, as required by rotational invari- 
ance. It is regularized in both methods by using two 
PV-type subtractions, involving both the fermion and 
boson propagators. 

This close connection between TLRS and PV-type 
regularization is of particular interest in further non- 
perturbative calculations in LFD In the TLRS 
scheme indeed, there is no need for additional non 
physical components in the state vector describing any 
physical system. Moreover, there is no large mass 
scale limit to perform numerically. This may ren- 
der possible large scale calculations of nonperturba- 
tive relativistic bound state systems in LFD. 
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